home *** CD-ROM | disk | FTP | other *** search
/ IRIX Base Documentation 2002 November / SGI IRIX Base Documentation 2002 November.iso / usr / share / catman / p_man / cat3 / SCSL / ztrsyl.z / ztrsyl
Encoding:
Text File  |  2002-10-03  |  4.4 KB  |  133 lines
  1.  
  2.  
  3.  
  4. ZZZZTTTTRRRRSSSSYYYYLLLL((((3333SSSS))))                                                          ZZZZTTTTRRRRSSSSYYYYLLLL((((3333SSSS))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      ZTRSYL - solve the complex Sylvester matrix equation
  10.  
  11. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  12.      SUBROUTINE ZTRSYL( TRANA, TRANB, ISGN, M, N, A, LDA, B, LDB, C, LDC,
  13.                         SCALE, INFO )
  14.  
  15.          CHARACTER      TRANA, TRANB
  16.  
  17.          INTEGER        INFO, ISGN, LDA, LDB, LDC, M, N
  18.  
  19.          DOUBLE         PRECISION SCALE
  20.  
  21.          COMPLEX*16     A( LDA, * ), B( LDB, * ), C( LDC, * )
  22.  
  23. IIIIMMMMPPPPLLLLEEEEMMMMEEEENNNNTTTTAAAATTTTIIIIOOOONNNN
  24.      These routines are part of the SCSL Scientific Library and can be loaded
  25.      using either the -lscs or the -lscs_mp option.  The -lscs_mp option
  26.      directs the linker to use the multi-processor version of the library.
  27.  
  28.      When linking to SCSL with -lscs or -lscs_mp, the default integer size is
  29.      4 bytes (32 bits). Another version of SCSL is available in which integers
  30.      are 8 bytes (64 bits).  This version allows the user access to larger
  31.      memory sizes and helps when porting legacy Cray codes.  It can be loaded
  32.      by using the -lscs_i8 option or the -lscs_i8_mp option. A program may use
  33.      only one of the two versions; 4-byte integer and 8-byte integer library
  34.      calls cannot be mixed.
  35.  
  36. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  37.      ZTRSYL solves the complex Sylvester matrix equation:
  38.         op(A)*X + X*op(B) = scale*C or
  39.         op(A)*X - X*op(B) = scale*C,
  40.  
  41.      where op(A) = A or A**H, and A and B are both upper triangular. A is M-
  42.      by-M and B is N-by-N; the right hand side C and the solution X are M-by-
  43.      N; and scale is an output scale factor, set <= 1 to avoid overflow in X.
  44.  
  45.  
  46. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  47.      TRANA   (input) CHARACTER*1
  48.              Specifies the option op(A):
  49.              = 'N': op(A) = A    (No transpose)
  50.              = 'C': op(A) = A**H (Conjugate transpose)
  51.  
  52.      TRANB   (input) CHARACTER*1
  53.              Specifies the option op(B):
  54.              = 'N': op(B) = B    (No transpose)
  55.              = 'C': op(B) = B**H (Conjugate transpose)
  56.  
  57.      ISGN    (input) INTEGER
  58.              Specifies the sign in the equation:
  59.              = +1: solve op(A)*X + X*op(B) = scale*C
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. ZZZZTTTTRRRRSSSSYYYYLLLL((((3333SSSS))))                                                          ZZZZTTTTRRRRSSSSYYYYLLLL((((3333SSSS))))
  71.  
  72.  
  73.  
  74.              = -1: solve op(A)*X - X*op(B) = scale*C
  75.  
  76.      M       (input) INTEGER
  77.              The order of the matrix A, and the number of rows in the matrices
  78.              X and C. M >= 0.
  79.  
  80.      N       (input) INTEGER
  81.              The order of the matrix B, and the number of columns in the
  82.              matrices X and C. N >= 0.
  83.  
  84.      A       (input) COMPLEX*16 array, dimension (LDA,M)
  85.              The upper triangular matrix A.
  86.  
  87.      LDA     (input) INTEGER
  88.              The leading dimension of the array A. LDA >= max(1,M).
  89.  
  90.      B       (input) COMPLEX*16 array, dimension (LDB,N)
  91.              The upper triangular matrix B.
  92.  
  93.      LDB     (input) INTEGER
  94.              The leading dimension of the array B. LDB >= max(1,N).
  95.  
  96.      C       (input/output) COMPLEX*16 array, dimension (LDC,N)
  97.              On entry, the M-by-N right hand side matrix C.  On exit, C is
  98.              overwritten by the solution matrix X.
  99.  
  100.      LDC     (input) INTEGER
  101.              The leading dimension of the array C. LDC >= max(1,M)
  102.  
  103.      SCALE   (output) DOUBLE PRECISION
  104.              The scale factor, scale, set <= 1 to avoid overflow in X.
  105.  
  106.      INFO    (output) INTEGER
  107.              = 0: successful exit
  108.              < 0: if INFO = -i, the i-th argument had an illegal value
  109.              = 1: A and B have common or very close eigenvalues; perturbed
  110.              values were used to solve the equation (but the matrices A and B
  111.              are unchanged).
  112.  
  113. SSSSEEEEEEEE AAAALLLLSSSSOOOO
  114.      INTRO_LAPACK(3S), INTRO_SCSL(3S)
  115.  
  116.      This man page is available only online.
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.